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ABSTRACT 


The  electrostatic  charge  on  airplanes  was  determined  from 
measurements  of  the  charge  induced  on  a  large  screen  while  the  air¬ 
planes  flew  over  it.  The  airplane  charge  was  computed  from  a  trans¬ 
fer  function  that  involved  the  relaxation  time  of  the  sensing  circuit 
and  the  passage  time  of  the  airplane.  Corrections  for  variation  in 
the  height  and  displacement  of  the  airplanes  were  made  from  labora¬ 
tory  measurements  on  scale  models  of  the  field  equipment.  The 
average  charge  on  a  jet-propelled  F-86  fighter  was  found  to  be 
-35.4  x  10-6  coulomb  and  on  a  B-26  bomber ,  -0.56  x  10~8  coulomb. 

The  corresponding  aircraft  potentials  are  approximately  -141,000  and 
-1200  v. 

1.  INTRODUCTION 

It  has  long  been  known  that  aircraft  in  flight  carry  an 
electric  charge.  The  effects  of  this  charge  are  usually  observed 
when  it  is  so  large  as  to  be  hazardous  as  when  the  Zeppelin 
Hindenburg  burst  into  flame  or  so  troublesome  as  when  precipita¬ 
tion  static  interferes  with  aircraft  communication  (ref  1,2).  Such 
phenomena  occur  when  the  static  charge  has  built  up  to  the  level 
where  a  spark  or  corona  occurs.  The  extent  of  charging  when  con¬ 
ditions  are  not  conducive  to  such  a  large  accumulation  is  not  so 
easily  recognized. 

The  present  investigation  was  made  to  determine  to  what  extent 
airplanes  are  charged  in  clear  weather. 

2 .  TEST  METHOD 


There  are  many  possible  methods  for  measuring  the  electrostatic 
charge  on  aircraft.  The  method  used  by  Ross  Gunn  and  his  associates 
(ref  1),  wherein  electric  field  meters  were  installed  at  various 
points  on  the  surface  of  an  airplane,  was  not  feasible  in  this  investi¬ 
gation  because  of  the  required  major  modifications  '„o  the  airplane. 
Flying  an  airplane  through  a  Faraday  cage  or  into  a  screen 
in  order  to  measure  the  full  charge  by  induction  or  by  charge  transfer 
was  obviously  not  feasible.  The  method  selected,  a  modification  of 
the  Faraday  cage  technique,  consisted  in  flying  the  airplanes  over 
a  large  plate  which  detected  the  charge  by  induction.  The  flights 
were  therefore  restricted  to  low  altitudes. 

The  test  equipment  was  arranged  as  shown  schematically  in 
figure  1.  P  represents  the  induction  plate  and  E  the  earth.  The  net 
resistance  and  capacitance  between  P  and  E  are  represented  by  R^  +1^  and 
Cg ,  respectively.  The  electric  charge  on  an  airplane,  A,  flying 
overhead  induces  a  charge  on  P.  A  voltage  proportional  to  that  on 
P  is  recorded  by  an  oscillograph,  B. 
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A  relationship  between  the  charge  and  the  recorded  signal  is 
derived  in  appendix  A.  The  equation  used  is 


in  which 


Q  = 


2C.  v 
2  max 

k  f(t) 


(1) 


max 


Q  =  charge  on  aircraft,  coul 

C  =  capacitance  between  plate  and  ground,  0.0236  pf 


k  =  proportionality  factor 
v  =  recorded  signal,  volts 
f(t)  =  transfer  function 

The  determination  of  the  proportionality  factor,  which  represents  the 
ratio  of  the  peak  induced  charge  on  the  plate  to  Q,  is  given  in 
appendix  B. 


3.  TEST  EQUIPMENT 


The  detection  equipment  (fig.  2)  was  set  up  on  an  apron  at  Phillips 
Field,  Aberdeen  Proving  Ground,  Maryland.  A  metal  ground  screen  con¬ 
sisting  of  2-in.  mesh  netting,  150  by  150  ft,  was  spread  out  to 
establish  a  definite  electrical  ground.  The  induction  plate  consisted 
of  50  wood  frames  covered  with  1-in.  mesh  wire  netting  and  supported 
by  1-ft  plastic  legs.  The  frames  were  assembled  over  the  center  of 
the  ground  screen  and  connected  together  to  form  a  plate  80  ft  wide 
by  55  ft  in  the  direction  of  the  aircraft  flight. 

The  recording  equipment  was  located  about  400  ft  to  the  side  of 
the  detection  equipment  (fig.  3) .  A  shielded  cable  connected  the  in¬ 
duction  plate  to  the  recorder.  The  capacitance  of  the  plate  and 
cable  to  ground  (Gg,  fig. 1)  was  measured  and  found  to  be  0.0236  pf. 

The  resistance  consisted  of  the  recorder  resistance  (Rg,  10  meg)  plus 
additional  series  resistance  (f^  )  needed  for  attenuation. 

Auxiliary  equipment  was  used  to  measure  the  velocity,  height, 
and  off-center  displacement  of  an  airplane  each  time  it  passed  over 
the  screen.  Two  photocell  instruments  were  placed  along  the  line  of 
flight,  400  ft  apart,  one  on  each  side  of  the  screen.  The  aircraft 
speed  was  computed  from  the  elapsed  time  between  the  two  signals 
caused  by  obscuration.  An  open  shutter  movie  camera  was  placed 
slightly  beyond  the  second  photocell  to  photograph  the  airplanes. 

The  height  and  position  of  each  airplane  was  determined  accurately 
from  the  resulting  photograph,  as  explained  in  appendix  C.  An  alidade 
located  to  the  side  of  the  screens  enabled  a  ground  observer  to 
give  the  pilot  immediate  information  on  the  height  of  passage  of  the 
airplane . 
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FIG.  2  PILOT'S  VIEW  OF  INDUCTION  PLATES 
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4. 


TEST  PROCEDURE 


The  electric  charge  was  measured  in  clear  weather  on  B-26  pro¬ 
peller-driven  bombers  flown  by  crews  of  the  6570th  Test  Group  of 
the  Air  Force  and  on  F-86  jet-propelled  figiters  flown  by  members 
of  the  Maryland  Air  National  Guard.  In  addition,  an  L-19  observa¬ 
tion  plane  and  a  C-54  cargo  plane  were  invited  to  fly  over  the 
induction  screen  for  observation  while  the  apparatus  was  in  place. 

The  pilots  flew  over  the  screen  as  low  as  safety  limits  per¬ 
mitted,  levelling  off  about  4000  ft  before  the  test  area  and,  guided 
by  ground  markers  and  by  the  edge  of  the  runway,  tried  to  fly  over 
the  center  of  the  screen.  Heights  over  the  screen  ranged  from  20  to 
63  ft,  and  displacements  from  the  center  line  ranged  from  zero  to 
25  ft.  The  velocities  were  usually  close  to  their  upper  limits  for 
the  combat  aircraft.  Original  plans  to  cover  a  range  of  velocities 
and  other  flight  situations  had  to  be  abandoned  because  the  permis¬ 
sible  flying  time  of  the  aircraft  was  cut  short. 

Recordings  of  the  induced  voltages,  such  as  those  shown  in 
figure  4,  were  obtained  for  each  airplane  passage.  The  upper  trace 
(fig.  4A)  is  typical  of  those  originally  obtained  for  a  B-26  with 
the  test  circuit  shown  in  figure  1  with  R.  =0.  The  very  good 
agreement  of  this  curve  with  the  theoretical  curve  A  in  figure  A3 
confirms  the  derivation  of  the  transfer  function  in  appendix  A. 
(Polarities  are  inverted.)  Inasmuch  as  the  only  data  desired  for 
equation  1  were  the  peak  values,  a  vacuum  tube  diode  was  later 
inserted  in  the  circuit  at  point  X  (fig.  1)  to  eliminate  the  signal 
reversal  and  so  establish  a  better  base  line.  The  resulting  trace 
for  a  B-26  is  shown  in  figure  4B  and  for  an  F-86  in  figure  4C. 

These  curves  are  for  conditions  comparable  to  those  of  curves  B  and 
C  in  figure  A3. 

5 .  RESULTS 


The  electrostatic  charges  on  the  airplanes  were  negative.  The 

magnitudes  were  computed  from  equation  1  for  those  passages  in 

which  complete  data  were  obtained.  The  heights  and  displacements 

were  obtained  in  accordance  with  the  procedure  outlined  in  appendix 

C,and  the  k  factor  was  obtained  in  accordance  with  the  procedure 

outlined  in  appendix  B.  The  peak  value  of  the  transfer  function 

f(t)  was  obtained  from  appendix  A.  The  peak  voltages  (v  ) 
max  b  max' 

were  read  from  the  oscillograph  traces.  The  results  are  tabulated 

for  the  F-86  and  B-26  in  tables  1  and  2,  respectively. 

The  average  electric  charge  found  on  the  F-86  was  -35.4  x  10-6 
coul  and  on  the  B-26,  -0.56  x  lO^coul.  The  signals  made  by  the 
passages  of  the  C-54  and  L-19  aircraft  were  similar  to  those  shown 
in  figure  4, and  their  average  amplitudes  were  -5.8  and  -4.6  v, 
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A.  B-26  AIRPLANE,  R,  =  0 


B.  B-26  AIRPLANE,  R,  =  100  MEGOHMS 


C.  F-86  AIRPLANE,  R,  =  100  MEGOHMS 

FIG.  4  RECORDINGS  OF  INDUCED  VOLTAGE 
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respectively c  Although  these  signals  exceeded  those  obtained  with 
the  B-26  bomber,  the  magnitudes  of  the  associated  airplane  charge 
could  not  be  evaluated  because  the  auxiliary  measurements  for 
height  and  displacement  were  not  recorded  properly. 

The  potential  of  an  airplane  may  be  computed  from  its  electric 
charge  if  its  capacitance  is  known.  Ross  Gunn  (ref  1)  gave  the 
capacitance  of  an  airplane  in  esu  as  being  approximately  equal  to 
20  percent  of  its  wing  span  in  cm.  In  accordance  with  that  criterion, 
the  capacitance  of  the  F-86  is 

(20$)  (37.1  ft)  (30.48  cm/ft)  (1.11  pf/cm)  =  251  pf 
The  average  potential  was  therefore 

(-35 . 4xl0_s  coul)/(251xl0"12  f)  =  -141,000  v. 

The  corresponding  values  for  the  B-26  are 

(20$)  (70  ft)  (30.48  cm/ft)  (1.11  pf/cm)  =  473  pf 


and 

(-0.56x10-6  coul)/(473xl0"12  f)  =  1200  v. 

The  magnitudes  of  the  charge  measured  on  the  aircraft  are 
shown  in  tables  1  and  2.  Of  20  readings  on  the  F-86,  only  one  ex¬ 
ceeded  52.4  pcoul,  and  only  one  was  less  than  26.3  pcoul.  For  the 
B-26,  four  readings  exceeded  0.85  pcoul,  five  were  less  than  0.27 
pcoul,  and  the  remaining  16  were  between  those  values.  These  dis¬ 
tributions  may  be  attributed  to  variation  in  the  electrostatic 
charge  on  the  aircraft  from  one  run  to  another  and  from  day  to  day, 

and  to  departure  of  the  true  values  of  k  and  f(t)  from  the  values 

max 

used  in  the  computations. 

Most  of  the  runs  were  made  at  maximum  flight  speeds.  The 
number  of  flights  at  lower  speeds  was  insufficient  to  determine 
any  significant  difference  in  the  charging  effects.  Any  possible 
effects  may  well  have  been  hidden  by  other  factors,  such  as  the  length 
of  time  flown  at  the  lower  speeds  or  passage  through  clouds  before 
passage  over  the  induction  screen.  No  significant  difference  was 
found  in  the  charge  on  the  B-26  between  runs  with  wicks  and  those 
without  wicks.  However,  radio  communications  between  the  airplane 
and  the  ground  were  nearly  severed  with  the  wicks  disconnected,  which 
indicated  that  the  airplane  was  losing  charge  in  a  very  deleterious 
manner,  perhaps  from  the  antenna  itself.  The  F-86  was  not  equipped 
with  wicks;  the  jet  exhaust  probably  served  the  purpose. 
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6. 


CONCLUSION 


The  electric  charge  on  various  military  aircraft  was  measured 
in  clear  weather  and  at  low  altitudes  by  a  method  of  charge  induc¬ 
tion  without  altering  the  construction  or  the  usual  flight  procedures 
of  the  aircraft.  The  polarity  of  the  charge  was  found  to  be  negative. 
The  average  charge  on  the  F-86  jet  aircraft,  35.4  pcoui,  was  much 
greater  than  that  on  the  B-26  propeller  aircraft,  0.56  pcoul.  These 
corresponded  to  potentials  on  the  aircraft  of  about  141,000  and 
1200  v,  respectively.  The  magnitudes  of  the  electric  charge  on 
L-19  and  C-54  airplanes  could  not  be  computed  although  they  were 
judged  to  be  significant  from  the  induced  voltage  measurements.  It 
would  be  desirable  to  conduct  further  tests  to  study  the  accumulation 
of  electric  charge  on  aircraft  under  different  conditions. 
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TABLE  1.  RESULTS  WITH  F-86  AIRPLANES* 


Velocity 

(mph) 

Height 

(#?S)b 

Displace¬ 

ment 

(#?S)b 

Proportion¬ 

ality 

(k) 

(vr2)c 

(meg) 

Peak 

Signal 

(v  v) 

max, 

Charge  on 
Airplane'* 
(10-«coul) 

590 

149 

32 

0.17 

110 

220 

32.0 

590 

127 

5 

0.24 

110 

275 

28.4 

345 

74 

38 

0.39 

110 

140 

8.9 

345 

66 

30 

0.44 

110 

465 

26.1 

580 

90 

46 

0.32 

110 

385 

29.8 

580 

138 

32 

0.20 

110 

275 

34.0 

580 

116 

5 

0.27 

110 

330 

30.2 

580 

119 

14 

0.26 

110 

360 

34.3 

580 

116 

24 

0.26 

110 

330 

31.4 

535 

100 

11 

0.32 

110 

440 

34.0 

560 

105 

19 

0.30 

110 

415 

34.2 

610 

54 

68 

0.43 

110 

910 

52.4 

600 

105 

23 

0.30 

110 

690 

56.9 

590 

72 

11 

0.42 

110 

495 

29.2 

610 

103 

50 

0.28 

110 

385 

34.0 

550 

130 

20 

0.22 

110 

410 

46.1 

560 

122 

24 

0.24 

110 

440 

45.4 

590 

66 

22 

0.44 

110 

700 

39.4 

610 

97 

22 

0.32 

110 

550 

42.5 

600 

77 

32 

0.38 

110 

605 

39.4 

Mean  35.4 


Notes : 

a.  Diode  in  circuit  for  all  runs. 

b.  The  height  of  the  lower  fuselage  surface  above  the  induction 
screen  and  the  displacement  of  the  flight  path  from  the  center  line  of 
the  plate  are  given  in  percentage  of  the  wing  span  (WS) .  The  wing  span 
of  the  F-86  is  37.1  ft. 

c.  R^+Rg  refers  to  the  circuit  of  figure  1  and  affects  the 

transfer  function.  f(t).  For  these  passages  f(t)  is  1.908; 

'  max 

d.  The  polarity  of  the  signal  voltage  and  the  charge  is 
negative . 


Nc 

ti 

li 

TI 

ti 

f< 
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TV  ^LE  2.  RESULTS  WITH  B-26  AIRPLANES 


Velocity 

(mph) 

Height 

($WS)a 

Displace 

ment 

(<*WS)a 

!-  Proportion-  i 
ality 

GO 

(VR2) 

(meg) 

b 

Peak 

Signal 

(vmax/> 

Charge  on  c 
Airplane 
(10-6  coul) 

330 

50 

28 

0.145 

10 

0.95 

0.50 

330 

90 

50 

0.065 

10 

1.65 

2.00 

320 

68 

45 

0.090 

10 

1.05 

0.85 

370 

65 

26 

0.110 

10 

0.75 

0.52 

340d 

5C 

11 

0.145 

10 

3.0 

1.59 

330 

57 

9 

0.110 

10 

2.6 

1.81 

320 

66 

14 

0.110 

10 

1.8 

1.28 

340 

68 

19 

0.105 

10 

0.6 

0.41 

340 

44 

32 

0.165 

10 

0.2 

0.09 

210 

32 

13 

0.245 

110 

3.3 

0.38 

215 

42 

6 

0.195 

110 

0.6 

0.09 

210 

41 

14 

0.195 

110 

3.6 

0.51 

260 

59 

7 

0.130 

110 

1.4 

0.30 

260 

34 

14 

0.230 

110 

2.2 

0.27 

315 

36 

19 

0.215 

110 

0.8 

0.11 

330e 

54 

11 

0.140 

110 

2.5 

0.49 

340 

53 

17 

0.140 

110 

1.7 

0.33 

340 

44 

9 

0.180 

110 

0.6 

0.09 

325 

33 

0 

0.250 

110 

2.8 

0.31 

310 

34 

0 

0.240 

110 

2.8 

0.32 

340 

44 

9 

0.180 

110 

0.3 

0.05 

325 

46 

10 

0.170 

110 

2.2 

0.33 

310 

46 

16 

0.170 

110 

2.2 

0.35 

300 

35 

0 

0.235 

110 

3.3 

0.40 

Mean 

0.56 

Notes : 

a. 

The  height  of  the  lower  fuselage 

i  surface 

above  the  indue- 

tion  screen  and  the  displacement  of  the  flight  path  from  the  center 

line  of  the  plate 

are  given  in  percentage 

of  the  wing  span 

(<£ws) . 

The  wing 

span  of  the  B-26 

airplane  is  70.0 

ft. 

b. 

R.  +  R_  refers  to 

the  circuit  of 

figure 

1 

and  affects  the 

transfer 

function, 

f(t). 

f(t)  =  0.616 

for  R, 

+ 

R-  =  10 

meg  and 

£<‘W 

1.677  for 

ni+Ra  = 

110  meg  for  these  passages. 

c. 

The  polarity  of 

the  signal  voltage  and 

the  charge 

is  negative. 

d. 

Diode  in 

circuit 

for  this  run  and  following  runs. 

e. 

Antistatic  wicks 

disconnected  for  this 

run  and  following  runs. 
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APPENDIX  A— DERIVATION  AND  EVALUATION  OF  THE  TRANSFER  FUNCTION 


The  induction  screen  over  which  the  airplane  flies  will  be  con¬ 
sidered  a  solid  plate,  P,  (fig.  1).  A  voltage  pulse  is  recorded  on 
the  oscillograph,  B,  for  each  flyover.  This  pulse  is  the  result  of 
an  electric  charge  q  induced  in  the  capacitor  and  the  capacitance 
Cj .  The  maximum  value  of  q  is 

qm  =  kQ  (Al) 

max 

where  Q  is  the  charge  on  the  airplane  and  0  <  k  <  1.  This  pro¬ 
portionality  factor,  k,  may  be  obtained  from  scale  model  measure¬ 
ments,  as  described  in  appendix  B. 

During  the  flyover,  the  charge  q  depends  on  an 
voltage  v^  (fig.  Al,  derived  from  fig.  1)  so  that 

q  =  -  v 

cl+ca  ° 

if  R  is  infinitely  large.  Also, 

q  =  qv  (A3) 

If  no  charge  were  lost  from  Cg  during  the  flyover,  q^  could  be 

computed  simply  from  v  read  on  the  oscillograph  tape.  However, 

since  charge  is  lost  through  R,  a  functional  relationship  between  q 
and  v  must  be  derived. 

The  characteristics  of  the  voltage  v  are  that  it  is  initially 
zero,  reaches  a  maximum  when  the  airplane  is  directly  over  the  screen, 
and  returns  symmetrically  to  zero  afterward,  as  indicated  in  figure  A2, 
It  may  be  considered  as  one  cycle  of  the  sinusoidal  wave  given  by 

(1  “  cos  0  <  t  <  T  (A4) 

Vo  "  |  0,  t  <  0,  t  >  T 

where  V  =  undetermined  maximum  value  of  v^ 
o  o 

t  =  elapsed  time 

T  =  period  during  which  the  charge  on  the  airplane  is  effective. 
From  figure  Al 


assumed  driving 

(A2) 


v 

o 


1 


C: 


i0dT  +  v 


(A5) 
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V  = 


j_»  ^ =  iB 


(A6) 


i  =  i  +  ia 
o 


(A7) 


From  equations  A4,  A5,  and  A6 


kU 


i-dj  +  iR  =  < 


0,  t  <  0 

(1  "  cos  0  <  t  <  T 


(A8) 


Zc 

2 


2nt> 

T 


r  'O  j  1  2TT(t-T)~l  , 

(l-  COS  — )  -  —  l^l-cos  - - -  j,  t>T 


Since  v  svsi  sisig  =  0  for  t  <  0,  the  application  of  Laplace 
transforms  yields 


^  ( - — — ),  0  <  t  <  T 

2s  ss  +  £2 


q 


+  R  i  =  > 


^  (- 
2s  v 


- - )  (1  ~e~TS),  t>T 


s2  +  P2 


in  which 


V 


p 


2n 


(A9) 


(A10) 


From  equations  A6  and  A7 } 


i  =  i  +  U 
o  * 


=  i(l  +  RQ, s) 


(All) 


Substituting  equation  All  into  A9  and  clearing  of  fractions 

I  vGq  e3  1 

r  "|  —  I  2  ^  sa  +  0<t  ^ 

[i  +  R<q  +  q>s] 1  =  ,  9  _TS 

L  J  V„  C,  B2  .]  -  e  Ts. 


V 


s2  + 


)j  t>T. 


(A12) 


In  the  interval  0<t<T 


-  Vo^2Y  !  ! 

i=— 2 - TT? 


(A13) 
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V 


i^dT  =  iR 


(A6) 


J_0 


i  =  i  +  ig 
o 

From  equations  A4,  A5,  and  A6 

(o,  t  <  0 
Vo  .. 


k  L* 


+  iR  = 


2nt 
T  J 

2nt. 


v  "TT  O 

-f  (1-  COS  —)  -  -j  (_ 


o 

A 

t  <  T 

V  r 

2  L 

1-cos 

2n(t-T)~ 

T 


Since  v  =v=i  =  i  =  i_  =  0  for  t  <  0,  the  application  of 
transforms  yields 

(y  os 

2§  (  a  0  <  t  <T 

2S  s3  +  p3 

P3 


+  R  i  =  ) 


in  which 


V 


2s  g2  +  ^2 

o  2rr 
P  = 


-)  (1  -e  TS),  t>T 


Fi'om  equations  A6  and  A7 , 

i  =  i  +  i, 
o  * 

=  i  ( 1  +  RCgS) 

Substituting  equation  All  into  A9  and  clearing  of  fractions 

p3 

i  - 

[i  +  R(q  +  Cg )  s']  T  = 


0<t  <T 


vc  qp*  ]  - 


-Ts 


V 


s3  +  p2 


-),  t>T. 


In  the  interval  0<t<T 


-  Vo9iP8Y 

i  =  - - - (- 


„  1  w  1 _ X 

s3  +  P2  n  S  +  y> 


(A7) 

(A8) 

L  t>T 
Laplace 

(A9) 

(A10) 

(All) 

(A12) 


(A13) 
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m  which 


_ 1 

Y  "  R(q  +  q) 


(A14) 


On  taking  the  inverse  Laplace  transform 


voqB3v  e~yt 

sin(Bt  -  (r) 

1  ~ 

2  B3  +  y8 

B(B3  +  y3^ 

in  which 

%  =  tan-1 

B 

V 

(  A15) 


(A16) 


In  the  interval  t>  T 


T  _  voC,  1-  e~TS 

2  (s3  +  3s)  (s  +  Y)l 


The  inverse  transform  is 


(A17) 


i  =  Wll  (_ 


8' 


B3  +  V3 


-)  (1  -  e^T)  e"^ 


(A18) 


Multiplying  equations  A15  or  A18  (according  to  whether  t<T  or 
t>T)  by  R  gives  the  measured  signal  v.  Rearranging  terms  gives 


v  = 


Cl 


2  q  +  Cg 


f(t) 


(A19) 


where 


f(t) 


r 


Yt+  (— air*8 ^  sin <Bt  -  t> 


B 


,  o<t<r 


b3  f/i  vT\  -  vtn 

— -  <1  -  eT  )  e  Y  , 


t>T 


(A20) 


B3  +  y2  L 


Equation  A19  includes  the  undetermined  constant  VQ  and  the  un¬ 
known  capacitance  q,  both  of  which  can  be  eliminated.  From  equation 
A2 
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Vo 


(  A21) 


max  q  +Cg 

Substituting  Vo  from  this  equation  into  Equation  A19 


v  = 


q  f  (t) 

max 


Substituting  from  equation  A1 


kQ  f(t) 

20, 


(  A22) 


(A23) 


from  which 


2Ca  v 
k  f(t) 


(  A24) 


Since  v  and  f(t)  attain  their  maximum  values  simultaneously 


Q  = 


2C_v 
2  max 

k  f(t) 


max 


(1) 


This  is  the  formula  given  in  the  text  for  the  computation  of  Q.  It 
can  be  seen  that  it  has  the  form  of  qv  multiplied  by  correction 
terms . 

Values  of  the  transfer  function  were  computed  not  only  at  its 
maximum  but  over  a  large  interval  of  time  in  order  to  check  its 
validity.  The  parameters,  defined  above  are 


0  = 


2tt 

T 


(A10) 


_ 1 

Y  "  R(q  +q> 


(  A14) 


i| i  =  tan  *  — —  (A16) 


T  =  effective  flyover  interval  (A4) 

R  =  10  meg  for  some  flights  and  110  meg  for  others,  q  was  measured 
and  found  to  be  0.0236  pf.  q  was  not  measured,  but  it  is  very 
much  smaller  than  q .  The  capacitance  of  the  B-26  bomber  to  all 


in  which 


_ 1 

Y  =  R(q  +  Cg) 

On  taking  the  inverse  Laplace  transform 

e”Yt  ~ 

_  sin(gt  -jr) 

P3  +  /  +  0(&3  +  ^)| 

if  =  tan-1  — 

Y 

In  the  interval  t>  T 


1  = 


in  which 


(A14) 


(A15) 


(A16) 


-Ts 


t  -  vq^  e2v  Iz _ _ 

2  (s3  +  p3)  (S  +  y)\ 


The  inverse  transform  is 


(A17) 


V^e2  v  .  e: 
1  = -  c - 


32  +  y2 


-)  (1  -  e^)  e"** 


(  A18) 


Multiplying  equations  A15  or  A18  (according  to  whether  t<T  or 
t>T)  by  R  gives  the  measured  signal  v.  Rearranging  terms  gives 


V  -- 


2  Ci  +  Ca 


-  f(t) 


(A19) 


where 


f(t) 


_  3a  _ 

32  +  / 

38 


;  Yt+  (S—i-Y2)!  sin(3t  - 
f(l  -  eYT)  e  ”  Y1  ,  ;  t>T 


,  0<t<T 


(A20) 


e3  +  y3  l 

Equation  A19  includes  the  undetermined  constant  VQ  and  the  un¬ 
known  capacitance  Cj ,  both  of  which  can  be  eliminated.  From  equation 
A2 
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(  A21) 


Vo 


max  Cj  -fCg 

Substituting  Vo  from  this  equation  into  Equation  A19 


v  = 


q  a  f^t) 

max 


2^ 


(A22) 


Substituting  from  equation  A1 


v  = 


kQ  f(t) 
2ca 


(A23) 


from  which 


Q  = 


2ca 


k  f(t) 


(A24) 


Since  v  and  f(t)  attain  their  maximum  values  simultaneously 


Q  = 


2C_v 
2  max 

k  f(t) 


max 


(1) 


This  is  the  formula  given  in  the  text  for  the  computation  of  Q.  It 
can  be  seen  that  it  has  the  form  of  CgV  multiplied  by  correction 
terms . 

Values  of  the  transfer  function  were  computed  not  only  at  its 
maximum  but  over  a  large  interval  of  time  in  order  to  check  its 
validity.  The  parameters,  defined  above  are 


8  = 


2tt 

T 


(A10) 


1 

Y  ~  RCq  +Cg ) 


(A14) 


*  =  tan  *  — —  (A16) 

Y 

T  =  effective  flyover  interval  (A4) 

R  =  10  meg  for  some  flights  and  110  meg  for  others.  was  measured 
and  found  to  be  0.0236  pf.  was  not  measured,  but  it  is  very 
much  smaller  than  Cg .  The  capacitance  of  the  B-26  bomber  to  all 
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space  was  computed  to  be  475  pf  and  its  parallel  plate  capacitance  to 
the  induction  plate  when  directly  over  it  was  computed  to  be  82  pf. 
The  values  for  the  F-86  fighter  are  smaller,  250  pf  and  32  pf . 
respectively.  Hence,  Cl+Cs  was  taken  as  0.0237  pf  for  the  bomber  and 
0.0236  pf  for  the  fighter.  T  was  taken  as  2  sec  for  the  bomber, 
which  corresponds  to  440  ft  on  either  side  of  the  induction  plate 
at  a  speed  of  300  mph.  T  was  taken  as  1/2  sec  for  the  smaller  jet 
fighter,  which  corresponds  to  220  ft  on  either  side  of  the  induction 
plate  at  a  speed  of  600  mph. 


The  results  of  the  computations  are  shown  graphically  in  figure  A3. 
The  similarity  of  curve  A  of  figure  A3  with  the  upper  trace  of  figure 
4  is  apparent.  Likewise,  the  curve  C  of  figure  A3  is  similar  to  the 
lower  trace  of  figure  4,  although  a  discharge  diode  was  present  in 
the  field  test  circuit.  It  can  be  noted  that  the  area  under  the 
positive  portion  of  each  computed  curve  equals  the  area  under  the 
negative  portion,  which  is  to  be  expected  inasmuch  as  the  induction 
plate  returns  to  zero  potent: al  after  the  airplane  departs. 


The  maximum  possible  value  of  f(t)  occurs  when  R=«  which  makes 

v=0,  f(t)=l  -  cos  St,  an  f(t)  =2.  f(t)  then  becomes  the  same  as  the 
1  max 

driving  function,  vQ  (fig.  A3)  and  equation  1  reduces  to  the  ideal 
case 


Q  = 


C2V 


max 


k 


(A25) 


The  following  maximum  values  of  f(t),  used  in  the  computations 


for  Q,  were  calculated  from 

equations  A10,  A14, 

A16,  and  A20. 

Aircraft 

R(meg) 

T(sec) 

(^max 

B-26 

10 

2 

0.616 

B-26 

110 

2 

1.677 

F-86 

110 

1/2 

1.908 

— 

» 

- 

2.000 

It  can  be  seen 

that  the  theoretical  limit  value 

of  f(t)  is  almost 

reached  with  R 
jet  airplane. 

=110  megohms 

,  especially  in  the 

tests  with  the  F- 
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FIGURE  A3.  DRIVING  AND  TRANSFER  FUNCTIONS 


APPENDIX  B— DETERMINATION  OF  THE  PROPORTIONALITY  FACTOR  k 


The  proportionality  factor  k  is  the  largest  fraction  of  the 
total  aircraft  charge,  Q,  which  is  induced  on  the  induction  plate 
for  a  given  flight  path  of  the  aircraft.  Thus,  if  q  is  the  charge 
on  the  plate 


k  = 


Suax 

Q 


(Bl) 


The  proportionality  factor  was  therefore  determined  for  those  posi¬ 
tions  for  which  the  airplane  is  at  the  point  of  closest  approach 
to  the  induction  plate  for  a  given  flight  path,  so  that  q  is 
maximum. 


Values  of  k  were  measured  in  the  laboratory  on  scale  models  of 
the  equipment  used  in  the  field.  A  model  of  the  B-26  airplane  having 
scale  factor  of  67:1  was  assembled  and  coated  with  silver  paint. 

(The  wing  span  of  a  B-26  is  70  ft;  that  of  the  model  was  12  1/2  in.) 

An  induction  plate  and  a  ground  screen  were  built  to  the  same  scale 
and  placed  in  the  proper  positions,  as  shown  in  figure  Bl.  The  model 
was  suspended  above  the  screen  so  that  its  belly  could  be  placed  at 
known  heights  and  at  known  displacements  off  the  center  line.  (The 
ruler  in  the  photograph  was  not  present  during  test  measurements.) 

A  similar  mockup  was  made  for  the  F-86  model,  which  had  a  scale  of 
48:1.  (The  wing  span  of  an  F-86  is  35.1  ft;  that  of  the  model  9.25  in.) 

The  ground  screen  was  connected  electrically  to  earth.  An 
electrostatic  voltmeter  was  connected  between  the  ground  screen  and 
the  induction  plate  to  measure  the  electric  charge  induced  on  the 
plate.  The  capacitance  of  the  induction  plate  configuration  was 
measured  over  the  complete  deflection  range  of  the  voltmeter  so  that 
the  charge  could  be  calculated  from  the  voltage  reading. 


A  measurement  of  the  proportionality  factor  was  started  with  the 
induction  plate  at  zero  potential  and  the  airplane  at  a  distance  of 
more  than  5  ft  from  the  screens.  The  airplane  was  charged  to  a 
potential  of  approximately  5000  v  and  brought  to  a  predetermined 
position  over  the  screens.  The  associated  electrostatic  voltmeter 
reading  was  noted,  from  which  the  corresponding  value  of  induced 
charge  was  later  computed  from  the  formula  Q  =  CV.  The  full  charge 
on  the  airplane  was  then  measured  by  bringing  it  into  contact  with 
the  induction  plate.  The  ratio  of  the  two  charge  measurements 
yielded  the  proportionality  factor  k. 


The  values  of  k  determined  in  this  way  for  different  combina¬ 
tions  of  height  and  displacement  are  listed  in  tables  Bl  and  B2,  in 
which  the  position  coordinates  are  given  in  terms  of  airplane  wing 
spans.  This  table  was  used  in  the  computations  of  the  charge  on  both 
the  F-86  and  the  B-26  aircraft.  Interpolation  between  tabulated 
values  was  done  graphically. 
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iGURE  Bl.  TEST  ARRANGEMENT  FOR  MEASUREMENT  OF  PROPORTIONALITY  FACTOR. 


TABLE  Bl.  PROPORTIONALITY  FACTOR,  k,  FOR  DIFFERENT  POSITIONS  OF 
THE  F86  AIRCRAFT 


Height* 

Displacement 

3$; 

from  Center  Line 

1.00 

0.00 

0.25 

0.50 

0.75 

0.00 

1.000 

1.000 

1.000 

1.000 

1.000 

0.25 

0.685 

0.703 

0.678 

0 . 593 

0.461 

0.50 

0.554 

0.519 

0.490 

0.404 

0.357 

0.75 

0.400 

0.403 

0.352 

0.320 

0.269 

1.00 

0.321 

0.310 

0.294 

0.269 

0.224 

1.25 

0.250 

0.250 

0.218 

0.187 

0.177 

1.50 

0.191 

0.175 

0.152 

0.124 

0.133 

2.00 

0.111 

0.099 

0.109 

0.072 

0.073 

TABLE  B2.  PROPORTIONALITY  FACTOR,  k,  FOR  DIFFERENT  POSITIONS  OF 
THE  B26  AIRCRAFT 


Height* 

Displacement 

* 

from  Center  Line 

0.00 

0.25 

0.50 

0.75 

0.00 

1.000 

1.000 

1.000 

1.000 

0.25 

0.317 

0.314 

0.267 

0.116 

0.50 

0.161 

0.152 

0.116 

0.076 

0.75 

0.098 

0.087 

0.082 

0.052 

1.00 

0.066 

0.077 

0.056 

0.038 

♦Expressed  as  fractions  of  a  wing  span 


APPENDIX  C— DETERMINATION  OF  POSITION  OF  AIRCRAFT 


The  position  of  the  aircraft  after  it  had  flown  past  the  screen 
was  determined  accurately  by  an  open  shutter  35-mm  movie  camera 
developed  for  another  project  (ref  3).  The  camera  was  located 
about  250  ft  from  the  center  of  the  screens  along  the  flight  path. 

It  was  directed  vertically  upward  with  its  lens  at  the  same  level  as 
the  induction  plate.  The  film  travelled  immediately  behind  a  trans¬ 
verse  slit  in  the  image  plane  in  a  direction  opposite  to  that  of  the 
airplane. 

The  camera  operated  as  a  camera  obscura.  Each  point  on  the 
airplane  overhead  was  imaged  on  the  film  as  it  passed  through  the 
imaginary  plane  determined  by  the  slit  and  the  lens.  An  image 
point  remained  on  the  film  as  the  film  travelled  backward  with  a 
speed  of  about  3  1/2  ft/sec  while  the  airplane  flew  forward.  The 
resulting  photograph  appears  to  be  a  distorted  shadowgraph  of  the 
airplane  (fig.  Cl). 


The  position  and  the  velocity  of  an  airplane  can  be  computed 
from  its  photograph  and  the  known  geometry  of  the  system  by  the  use 
of  similar  triangles.  Figure  C2  indicates  how  this  is  accomplished. 
Since  dimensions  transverse  to  the  flight  path  are  not  distorted  by 
the  motion, 


h 

f 


a 

~b 


(Cl) 


where  h  =  unknown  height  of  a  photographed  section  of  the  airplane 

f  =  focal  length  of  the  camera  lens,  1.89  in.  (0.1575  ft) 

a  =  distance  between  two  known  transverse  points  on  the 
airplane,  obtained  from  dimension  drawing  of  the 
airplane 

b  =  measured  distance  between  the  two  corresponding  points 
on  the  photograph 

For  example,  the  diameter  of  the  image  of  the  engine  housing  of 
the  B-26  on  the  film  of  figure  Cl  is  0.315  in.  The  actual  diameter 
of  the  housing  is  60  in.  Hence, 


h  =  (0.1575)  ( 


60 


0.315 


-)  =  30  ft 


The  height  h  is  then  corrected,  again  by  means  of  the 
dimension  drawings;  to  obtain  the  height  of  the  belly  of  the  air¬ 
plane.  The  distance  the  airplane  center  is  off  the  desired  line  of 
flight  is  easily  determined  by  noting  which  part  of  the  image  lies 
in  the  center  of  the  film.  The  heights  and  the  displacements 
measured  by  this  procedure  are  given  in  tables  1  and  2  in  terms  of 
the  airplane  wing  span. 


B  —  26 


F-  86 

FIG.  Cl  PHOTOGRAPHS  OF  AIRPLANE  OVERHEAD 


Once  the  height  has  been  found  from  the  transverse  dimensions, 
the  airplane  velocity  can  be  found  from  the  longitudinal  dimensions 
and  the  film  velocity.  This  was  not  done  in  the  present  case  inas¬ 
much  as  the  velocity  was  measured  more  directly  with  two  photocells 
on  the  ground. 
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